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ABSTRACT 

In  this  work  some  approaches  used  In  classical  stat.is- 
tical  mechanics  for  deriving  Integral  equations  for  the  pair 
distribution  function  have  been  generalized  to  quantum 
systems.   In  the  functional  differentiation  method, 
distributions  depending  on  space  and  temperature  arguments 
are  generated  from  the  grand  partition  function  in  external 
field.   These  distributions  are  diagonal  elements  of  Green's 
functions  and  represent  physical  quantities  when  temperature 
arguments  are  set  equal.   A  given  distribution  can  be 
expanded  in  the  external  field  by  means  of  a  functional 
Taylor  series.   An  integral  equation  results  when  the  external 
field  is  specialized  to  one  arising  from  a  number  of  fixed 
particles  via  their  pair  interaction  with  the  remaining 
particles.   In  this  case  a  distribution  in  external  field 
may  be  converted  to  a  higher  order  distribution  for  a 
uniform  system.   This  process  is  accomplished  by  using  a 
Peynman  path  integral  representation  for  all  quantities. 
Classical  and  quantum  statistics  are  treated.   In  the  former 
case  the  problem  of  bounds  on  distributions  is  also  investi- 
gated.  An  alternate  approach,  the  potential  ensemble  method, 
applies  to  Debye-Huckel  sequences.   Here  one,  e.g.,  replaces 


ill 


the  actual  density  (with  a  particle  fixed)  by  one  In  which 
a  pair  Interaction  of  particles  Is  replaced  by  an  external 
field,  but  an  ensemble  average  must  be  taken  over  all 
external  fields.   A  steepest  descent  calculation  enables  one 
to  find  the  pair  distribution  and  best  external  field  self- 
conslstently.   The  quantum  formulation  of  this  method 
requires  all  quantities  to  be  expressed  on  a  lattice  In 
occupation  representation,  in  which  a  lattice  site  is 
specified  by  position  in  space,  temperature  and  particle 
number.   Agreement  between  the  two  methods  described  here 
is  indicated. 


iv 


1.   Introduction 

Techniques  exist  for  deriving  Integral  equations  for 
the  pair  distribution  In  classical  statistical  mechanics  [1]. 
In  one  of  these,  one  first  functionally  differentiates  the 
grand  partition  function  with  respect  to  an  external 
potential  generating  coordinate  distributions.   Then  one 
performs  a  functional  Taylor  expansion  between  two  quantities 
dependent  on  the  external  potential,  as  density  v.  potential. 
One  completes  the  derivation  by  using  a  relation  which 
connects  lower  to  higher  distributions.   This  relation 
Involves  "fixing"  one  or  more  particles.   In  this  fashion 
Integral  equations  of  the  Debye-Hiickel  [2],  Percus-Yevlck  [3], 
Klrkwood-Salsburg  [h]    and  Mayer-Montroll  [5]  types  can  be 
derived.   In  addition,  rigorous  bounds  on  distributions  are 
obtainable  by  truncating  the  latter  two  sequences,  retaining 
the  remainder  [1],  [6]. 

Another  technique  [7],  at  least  for  Debye-Huckel  sequences, 
permits  us  to  replace  the  actual  density  (tied  to  a  fixed 
particle)  by  one  in  which  the  internal  potential  is  replaced 
by  an  external  one  and  an  ensemble  average  is  taken  over  all 
external  potentials.   A  steepest  descent  calculation  enables 
one  to  find  the  pair  distribution  and  "best"  external 
potential  self-consistently. 

This  work  generalizes  the  above  to  quantum  systems. 
After  a  brief  survey  of  classical  theory,  we  begin  in  Section  5 
by  introducing  the  distributions  appropriate  to  quantum 


systems.   These  fimctlons  have  space  and  temperature 
arguments  but  are  less  general  than  Green's  functions,  being 
only  diagonal  elements  of  them.   The  observable  distributions 
are  found  by  setting  temperature  arguments  equal.  An 
Integral  representation  of  the  grand  partition  function  Is 
crucial  to  the  generalization  of  the  Taylor  series  method, 
and  Is  carried  out  via  the  Feynman  path  Integral  [8]  In 
Section  ^.   In  later  sections  we  see  to  what  extent  the 
tools  developed  can  provide  us  with  Integral  equations  and 
bounds  for  the  pair  distribution.   The  "potential  ensemble" 
[7],  [9]  method  and  Its  generalization  Is  discussed  In  the 
last  three  sections. 


2.   Survey  of  Classical  Theory 

The  canonical  partition  function  for  a  classical  N 
particle  system  is 

(2.1)  Zj^  =J    dx^  ..  dxjj  dp3_  ..  dp^  exp  -  m^   , 

with 

2 

N   pf 

(2.2)  Hj^=X:2?r-  +  |=^(^i'^-) 

and  P  =  (kT)""^.   The  grand  partition  function  is  a  weighted 
average 


(2.3)  3=2 


oo     ^PiiN 


N=0 


NJ         ^N 


which  upon  doing  the  momentum  Integrations  becomes 

00   N  p  

(2-^)   ^  =  ZZ|t-  /  dx   •.  dx   exp  -  P  5~  <f'(x.  ,x  ) 


with 


=  eP^^(_iIi^)5/2 


The  sth  particle  distribution  is  defined  as 

(2.5)   ng(xj  ..  Xg)  = 

.N 


^"^  IZtTtW  /^1--^N  exp  -  p|-:<i)(xi,x.) 


•  6(x^  -  ^i)  ••  5(Xg  -  Xg) 


I 


and  is  gotten  from^  by  appending  an  external  potential 

N 
5   U(x. )  to  it  and  taking  functional  derivatives.   In  general, 


1=1 


6e       . . . 5e 


Modified  distribution,  whose  arguments  can  refer  to 
Identical  particles  are  given  by 


5:3.[U] 


(2.7)   n^(x,  ..  xJU)^[U]  =  6_pu(;^1!.;5-0mx^)  '    ^'^" 


(2.8)   n(x)  =  n(x),  n^{x^x^)    =  n^{x^x^)    +   n(x-j^)5(x^  -  Xg)  , 


n^Cx-j^XgX^)  =  n^tx^XgX^)  +  n2(x^X2)5(x2  -  x^) 

+  n2(x2X^)6(x^  -  x^)  +  n2(x,x^)5(x^  -  x^) 
+  n(x^)5(x2  -  x^)5(x^  -  x^),... 

Ursell  functions,  which  vanish  whenever  one  particle  Is 
uncorrelated  with  the  others,  are  given  by 


PU(X,  )    PU(xJ 


(2.9)   Pg(x^  ..  Xg|U)  e    ^  ... 


6^  log  ^[U] ^ 

-PU(x^)      -PU(Xg) 
6e       .  .  .  6e 


e.g. 


(2.10)   F(x)  =  n(x),  F(x^X2)  =  n2(x-j^X2)  -  n(x-L)n(x2)  , 


F(x^X2X,)  =  n^(x^X2X,)  -  n{7i-^)n^{x^x    )    -   n(x2)n2(x-,^x,) 


-  n(x,)n2(x^X2)  +  2n(x^ )n(x2)n(x^) , . . . 


Finally,  modified  Ursell  functions 


6"  log  ^[U] 


(2.11)   F^U,  ..  xJU)  =  s-0U(.°f.^'.5%(xJ  '  -8 


A        A   /\ 


A   A.         _  A  A   ^\ 


(2.12)   F^  =  n^,  ^2  ^  ^2  "  "^1*^1'  ^3   ^  ^3   ~   ^       "^2^1  "^  Sn-j^n^n-j^, . .  .  , 


as  In  (2.10). 


Consider  now  a  t-partlcle  distribution  with  s  other 
particles  x,  ...  x  held  fixed.   The  fixed  particles  provide 
an  external  field  to  the  remaining  ones  because  of  the 
pair  interaction.   Explicitly, 

s 

(2.13)  U,(y)  =  YZ  My,x.)  . 

^         j=l         '^ 

It  can  be  shown  that 

/o  n,M  f  iTT  ^     '"s+t(^r-^s-yr--yt^ 

(2.14)  n^iV:^    ...    y^\\J^)    =   n  (x,  .  .  .x  J 


and 

(2.15)  ^[U,]  =  z"^  ng(x-^  ...  Xg)^  •  exp  P^^  M^^Xj)  . 

As  a  simple  application  of  (2.11),  (2.8)  and  (2.14),  the 
latter  with  s  =  t  =  1,  expand  n,(y|Ui )  in  a  functional 
Taylor  series  against  Ui . 

With  U.(y)  now  just  <t)(y,x), 

^o(y^)  r   5nfv) 

(2.16)  n^(y|U^)  =  4^3^  =  n(y)  +  j  g^J^  U^o  ^^^^  "  ^^^^  +  "" 

^Q^y^")  r  2 > , , 

(2.17)  -^ =  n  -  p  /  dz(n2(yz)  +  n5(y  -  z)  -  n  )^{z   -   x)  . 

This  is  the  Debye-Huckel  equation  [2],  solvable  by  Fourier 
transformation.   While  (2.17)  is  valid  for  weak  potentials. 


*  Note;   (2.14)  follows  from  (2.15)  by  applying  (2.6). 


other  choices  for  functlonals  in  the  Taylor  series  are 
Indicated  If  the  potentials  should  be  singular.   For  example 
ne^  V.  e~   to  first  order  yields  the  beginning  of  a 
Klrkwood-Salsburg  type  hierarchy. 

We  next  consider  a  way  of  finding  bounds  on  distributions 

-PU.(x) 
For  non-negative  potentials,  co(x)  =  e    ^    -  1,  (the 

Mayer-f  function),  Is  <_  0,  and 

(2.18)     /"^s^-^l  *•  ^s)^(^i)  ••  a3(Xg)dx-,_  ..  dXg 


has  the  sign  of  (-1)^,  as  distributions  are  non-negative. 
Taking  advantage  of  this  fact  we  expand 

n^(yi  ...  ytlU^)^tU^]  exp  P^U^(y.) 

-PU,  s 

against  e   ^.   Here  we  choose  Ui(y)  ==  >   (|)(y,x.).   With 

the  help  of  (2.6),  (2.14)  and  (2.15)  the  full  expansion  with 
remainder  is: 


(2.19)   ng_^^(x^  ..  Xg,y-L  ...y^)z"^  exp  p[y~  H^jy^^) +o5~"  <t>(Xj_Xj) 

i      1   /- 

=  ^^!  /  ^^1  •••  "^V  '^(^l^  •••  "^(^m) 

m=0    ^ 

^t+m^^l  •••  y^  ^1  •••  ^m^  +  ^sti  • 


We  do  not  give  the  remainder  explicitly.   Sufficient  to  say 
that  since  It  contains  a  product  of  ^  +  1   oor-factors  Its 
sign  is  simply  (-1) 

With  t  =  0,  s  varying, (2.19)  Is  the  Mayer-Montroll 
sequence.   With  s  =  1,  t  varying.  It  Is  the  Klrkwood-Salsburg 
sequence. 

The  Mayer-Montroll  sequence  has  the  form 


(2.20)     n.  =5=-'  (-!>"  ^m%  *  (-1)'^'  Rs^  • 


with  non-negative  Asm  and  ^   p,    and  no  z  dependence  In  the 
Asm.   The  i  =  0  term  Is 


(2.21)  n^  =  z^  A^^  -  R    ;    (n  =  1) 

s       so    so  ^    ^  o    ' 


the  £  =  1   term  with  n-,  from  (2.21) 


(2.22)  n^  =  z^  A„^  -  z^A^, (zA,   -  R,  )  +  R  ,  . 

s       so      sl^   lo    lo'    sl 

Next  i  =  2,  with  n^  from  (2.22)  and  n^  from  (2.21)  etc. 
In  this  way  an  alternating  sequence  of  upper  and  lower 
bounds  is  generated: 

(2.23)  n  <  z^A 

^    -^'     s  —    so 

>  z^A^^  -  z^+l  A^tA, 

—  so        sl  lo 

<  z^A    -  z^"*"-"-  A  ,At   +  z^"^^(A  .A,, A,   +  A  _A^  ), 

—  so         sl  lo       ^  sl  11  lo    s2  2o" 


these  being  truncations  of  the  full  Mayer  fugaclty  expansion. 


3.   Quantum  Partition  Functions  and  Distributions 

The  partition  function  for  an  N  particle  system 
obeying  Maxwell -Bolt zmann  (classical)  statistics  is 

(5.1)  Zn  =^  j  ^""l  ••  "^""n  5i(^i  ••  ^n^  ®     ^i^^'l  ••  "^N^ 

(3.2)  =  tr  e   ^  , 


with  no  symmetry  restrictions  on  the  wave  functions.   The 
grand  partition  function  is  the  same  weighted  average  as 
(2.3): 


00  ^N  ^^  ^-PH, 


(3.3)  ^  =ZZiT  t 

'^        N=0  ^• 

but  with  z  =  e^  from  here  on.   Alternatively,  if  a  "thermal 
propagator"  is  defined, 

(3.^)      KjjCx-L    ..    Xjj   p',y^    ...    yj^  p"  )    = 

=  zi$r(yi  ••  ^N^Ji^^i  ••  ^N^  •  ^^p  -  ^N,i(P'  -  P"  ^ 

p'    >  p" 
Kj^   =   0,     P'    <   P" 
then    (3*3)    goes   over  into 
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z^    r 

(5.5)  ^  -  ^  NT  J    ^^1  •  •  ^^N  ^N^^l  •  •  ^N  P'^1  •  •  ^N  °)  • 

The  designation  "propagator"  carries  over  from  quantum 

mechanics  where  an  identical  function  (except  for  p  — >  -rr) 

n 

advances  a  wave  function  in  time.   (See  also  Section  4.) 
The  one  particle  propagator  is 


(3.6)   K(xP',yP"  )  =  2Z  i>liy)^^i^)    exp  -  E^O'  -  P"  )  . 


Evaluating  for  free  particles 


p  P 
(5.7)   K°(xP',yP"  )  =IZ^  exp[lk(x  -  y)  -  ^^   O'  -  P'«  )] 


with  Fourier  transform 


P  P 
(3.8)       K°(k,P'  -  P"  )  =  exp  -  ^^  O'  -  P"  )  . 


If  the  particles  are  bosons  (fermions)  drop  the  NI  and  sum 
over  distinct  symmetric  (antisymmetric)  states  in  {J>.1)  or 
(3.4). 

Immediate  developments  will  not  depend  on  the  nature 
of  the  statistics.   We  define  s-partlcle  space-temperature 
distributions 


(3.9)   np'^(x{p{^2^2  •••  Vn^^n+lCl  •.•;•••;..•  ^^P^) 


-1  V—   zN 


:  6(x^OJ)  -  x{)  ...  6(x-l(P^)  -  Xj!^) 


I     I 


where  the  arguments  are  distributed  consecutively  Into  j 
groups  set  off  by  j_  .   The  J  labels  particles  1  ..  j. 
Within  each  group  arguments  refer  to  the  same  particle. 
If  j  =  s  for  example,  so  that  each  group  has  one  particle, 
the  distribution  Is  fully  mutual.   At  the  other  extreme 
J  =  1  and  the  distribution  Is  fully  self.   Each  5-functlon 
Is  a  Helsenberg  operator 

(3.10)  6(x(P')  -  xM  =  expO'H^)6(x  -  x')  exp  -  P'Hj^  , 

and  the  :  :  orders  these  factors  so  that  ones  with  larger  P' 
are  placed  further  to  the  left.   In  particular 

(3.11)  n(x'P')  -^  -^  ^|y  N  tr(e        ^6(x3^-x')e    ^), 
the  density,  a  constant,  n,  for  uniform  systems. 


10 


I    I    T    t 


(3.12)  n^{x^^-^x^^^)    =^ 


-1 


r  -0-Pi)H. 


N 


5(xi  - 


y-\) 


e  5{x^  -  X2)e 


o;  >  p^) 


the  self  pair  distribution,  and 


II       II 


-1 


.N 


(3.13)      Hgtx-^P-^iXgPg)    =^'     ZZfcN{N   -   Dtrle 


-  -(P-Pi)H. 


N 


&{Xt   - 


^1) 


e  5(X2  -  Xgje 


0{  >  ^2^  , 

the  mutual  pair  distribution.   The  physical  pair  distribution 
is  (3' 13)  with  p^  =  Pp  .   At  equal  temperature  arguments 
(3'12)  becomes  singular:   n5(x-,  -  x  ). 

We  append  to  -^   an  external  potential  in  the  Heisenberg 


<:r- 


representation 


(3.14)   ^[U]  =XI|t  *^  ^^P^-P%^  :  exp  -J   y~  U(x^(P'),P')dP' 


We  can  by  this  purely  formal  device  generate  quantum  n  and 
F  .   Equations  (2-7)  and  (2.8)  generalize  to 


*  Note:   (3.14)  is  to  be  regarded  as  in  the  interaction 

representation.   All  the  U  dependence  is  in  the  exponential, 
Equation  (3. 10)  is  still  true,  i.e.  these  operators  have 
no  U  dependence. 


11 


(3.15)      nCx-^P^L    • 


^s^s^^  ~   -5U(x.p  ")!  .  .-6U(x    p. 


u  =  0 


(3.16)      n(l)    =n{l),      n^{l2)    =n2(l;2)    +112(12)    , 


while 
(3.17) 


n^(123)    =  n^(l;2;3)    +  n^(l;23)    +  n^(13;2)    +  n^(12;3) 

+  n^(123),    ...    , 


P    fx    6  X   B    )    -  6    log>-[U] I 

^s^^l^l    •••    ^s^s^    -   -&U(x^fi-L)...-5U(Xgf^g)    |u  =   0   ' 


(3.18) 


'S  ^  ^  ^  ^^ 

F,    =  n,  ,      F„   =  n^   -  nn   , 


F^   =  n^   -  2. ^p^  +  2nnn,    ...    , 


generalize    (2.11)    and   (2.12).      To  close   this   section  some 
free  particle   functions   are   evaluated. 

For  classical   statistics    (3-5)    becomes: 


N      r  N 


(3.19)     ;i"   -ZIftJ    dx^    ...    dxj^  ]T^    K   (x^Px^O) 


=  YZwr  U    dx   K"(xPxO)) 


i=l 
N 


0 


(3.20)  =  exp   zV(-i^)^/2    ^ 

■2TrPTi^ 


by   (3.6), 

(3.20M    n°  =f     ^  ^°^  >"     =  -^-^^^^^    ' 


12 


The  free  particle  self  pair  distribution  from  (3.12) 


■1     z%  .    -P^N-l 


(5.21)   n°  =^      (^ll^tre   ^-^) 


-(P-Pt)t,        ,  -o'-p')T         ,  -p't-l 
tr(e     ^   -^  5{x-^  -  x^)e    -^   ^   -^  5(x^  _  y^^)^      ^   -^), 

which,  inserting  complete  sets  of  states. 


(3.22)      =  z  ;^ 


-0-p{)T        X*  1   , 


I   I  » 


r      ^*       '  T-     -(P-Pi)E. 
(3.23)    =  z  /  dx  XZli  (^l)Il  (^)e     ■"   ^ 


t 


J    ^ 


I   I 


.  -  —  j_   —  J.  ^ 

1 


L_  J  j^(x)5^(x2)e 

i 

(3.24)  =  z  Tdx  K°(xP,x|p-|^)K°(x|p]^,X2P2)K°(x2P2^xO)  . 
Then,  using  (3.7)  or  (5.8)  : 

(5.25)  n°  (k,p^-p^)  =ne"  ^    '"l"^^!-   3 
An  easy  extension  of  (3*24)  is: 
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(3.26)   n°(x^p-^  ...  XgPg)  =  z  Tdx  K°{xPx-lP-l)K°(x-^P-lX2P2) 


..  K°(XgpgXO)  , 


free  full  self  distributions,   (p>p>p>...>p  >0) 

The  free  functions  for  quantum  statistics  are  postponed  to 
a  later  section. 


k .      The  Integral  Representation 

Refer  back  to  the  expression  for  the  one  particle 
thermal  propagator,  (3.6).   If  an  arbitrary  wave  function 
V'(x,P)  is  defined  by 


(4.1)     ^(xp)  =  ^~c,     <j),  (x) 


¥ 


•k  ^k' 


le 


-^\ 


we  see  that  K  has  the  property 


(4.2)  fU^^^)^  J        K(x2P2^^p^)^(x^P^)dx-L  . 
Also  from  (3.6) 

(4.3)  K(x^p^x^p-^)  =J     K.{ji^^^x^^^)K{x^^^x^^^)ax^ 
Consider  now  the  following  expression: 


IK 


2 

(4.4)  I  =(_i2^)5/2  /    exp[-  -^-i—^ e(l)(x3_)]^(xiP^)dXi  , 

27reh     ^  -co       2h 

with  e  a  small  quantity  of  units  Inverse  energy  (like  P-,). 

When  Xp  Is  appreciably  different  from  x-,  little  Is 
contributed  to  the  Integral.   We  therefore  replace  x,  -  Xp 
by  T],  with  the  expectation  that  the  main  contribution  to 
(4.4)  occurs  at  small  t^. 


^  /^   00  2 

(4.5)   I  -   (-^)^/^  /    dri  exp[-  -^  ^  "  ^M^p  +  n )  ] 
.2-n-eti^     ^-oo          2Ti^  ^       ^ 

•  ^(x^  +  n^^i)  • 


Dropping  the  x]   In  the  potential  and  expanding  that  term 
in  e,  ^  In  T) : 


p 
(4.6)   I  =  (_m   )V2  /  ^^  g^p[_  -^n-.](l  -  e^x  )) 

P 
2         Sx^ 


(4.7)    =  (1  -  eMx2))^(x2Pi)  +|^e^  , 


dropping  the  smallest  term.   But 


<-4fr|;l^^(-2))^(-2Pi)=--^^   ' 
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Is  Schrodlnger ' s  equation  with  the  replacement 

-^  ->  Pi  so 

S^(x^p. ) 

(4.8)  I  =  ^(xg,?^)  +  e  ^1  ^   =  ^(xg,?!  +  e)  , 

to  first  order.   Hence  by  (4.2)  we  Identify 

(4.9)  K(x_p^    +  e,x   p   )    =   (^IL^)5/2  exp[-  -J^  (x     -  x    )2-   e^{^    )] 

"^   ^  ^   ^  27reh  2li^e 

Now  by  (M  -  1)  repeated  application  of  (4.3),  we  arrive  at 
the  representation,  with  some  notation  changes 


(4.10)   K(xP',yP"  )  -  a'W    dx^   ...  dx^ 

^  -oo    ^1       ^M-l 


M        .p 
exp  -  >   e[-rp  ^^  +  M^a  )]  > 
k^   2h^  Pr      Pk 


e  — >  0,  where  the  position  of  the  particle  at  "time'  p^^ 

is  denoted  by  Xg  ,  Pj^  -  ^^_-^   =  e,  Pj^  -  P^  =  Me, 

^k 

Pk     Pk  ^ 
(^.11)   ^p^  =      e    ^   '  Pm  =  P''  Po  =  P"  ' 

A   -  I  ^  \ ^/^ 

^  -  I  iT>"^   '  -^R'  ~-^>  ^ftti  ~y  • 

27reh         P        P 


This  Is  the  Feynman  path  Integral  form  of  the  thermal  propagator. 
We  may  envision  a  particle  traveling  from  y  at  "time" 


P"  to  X  at  P'  along  a  path  marked  off  by  the  Xg  ;  the  path 

^k 
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Is  weighted  by  the  exponential  in  (4.10)  and  all  x^  are 
integrated  over,  except  for  the  end  points. 

Equation  (4.11)  has  the  N  particle  generalization 


(4.12)   K(x^  ..  Xj^  P',y^  ..  y^  P"  )  = 


=  A     /  dXi   ...  dx,, 

Pi       Pm-1 


M         r  N  „  ,       N 

exp   -  XI  e  ^-r2  ^i      +5-2ZHXi      -X         ) 


k=l 


^1  2^2   .p^  .    i^j     .p^   .p^_, 


The  integral  representation  of  the  Boltzmann  grand  partition 
function  now  follows  from  (3.5) 

N  -  - 

(^•15)   J^  =  I~^^^J   dx-L    ..    dx^   /   dx        ...    dx  exp   -    [      ] 

^  Pi  PM-1 

(4.14)  Pm  =  P   '      Pq  =  ° 


X^        -    X^        -    X^,      .  .  .  ,     X  -    X^        -    Xjj     . 

Op  op 


Thus  each  Boltzmann  particle  returns  to  its  starting  place 
after  "time"  p. 

In  quantum  statistics  we  must  include  paths  in  which 
the  final  configuration  is  some  permutation  of  the  original 
one,  as  the  particles  are  indistinguishable  [8]. 

If  (x-,,...,x„)  denotes  the  original  configuration. 
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p(x-,    ...    x>t)    will   denote   the   final   one   and 


N 


(t.i5)   i=^|T«™5Z(ii)^/  <3=^i  •••  '3='N 

^        N=0  p  >-' 


dx-,       .  .  .    dXjj  exp   -    [       ]    , 

P  Pi  Pm-1 


/  meaning  a  particular  permutation  described  above.   The  + 

^P 
sign  is  for  bosons,  -  for  fermions. 

•  To  obtain  distributions  add  an  external  potential 


M      N 

(4.16)  u  =  XZ  e  .^_  u(x.  ,\) 

k^i   1-1   ^Pj^ 

to  «i  and  employ  (3.15)  -  (3-l8).   We  can  then  identify, 
for  classical  statistics, 

(4.17)  4J)(xX  •••  Vn'^n+ld  •••;•..;•••  ^^s^ 


Pm-1 


exp  -  [   ]     •  5(x,   -X-,  )  ...  5(x,   -X 

1  n 


5(x^    -X'   )  ...  5(x    -xM  , 
cgj     n+1         Ja'   ^ 
n+1  s 


t 
) 


\ 
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adding  notation  ^j^  (+  1)^,  /   In  quantum  statistics 

P         ^P 

See  (3-9)  and  discussion  following  It. 
Some  examples: 


(4.18,   M.'.n^y'iZj^'^'"']^^,...^.^ 


dx-,   .  .  .  dx„    exp  -  [   ] 

\  \_^  U  =  0 


•  5(x     -  X')  , 


(4.19)   n^Cx'P'x"  P"  )   ("self") 


1  •  •  ^^N 


dx,   .  .  dx,.    exp  -  [   ] 


6(x    -  x')5(x,    -  X"  )  , 


(4.20)   n(x'P»;x"P")    ("mutual") 


-I  N 

,-1  '^c z     «MN 


J.       >__Xn=27T^   ^  dx^  ..  dx^^  dx    ..  dx 

•^^1      PM-1 


exp  -  [   ]      •  5(x,   -  x')5(x^   -X"  )  , 
U  =  0        pi         '=^pn 
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adding  the  usual  notation  for  quantum  statistics.   It  Is 

also  true  that,  for  all  quantities  expressed  as  path  integrals, 

(4.21)   n^(x^P,  ...  VslU)>fU]  =  ,5U(xSt)>— 5U(x  p  )  ' 


and 


(4.22)   P3(x,P,...  x^P^IU)  =  _60(x^pi[)^?!.5U(.3P3)  ' 


w 


ith  the  U  of  (4.l6) 


The  n  have  also  the  simple  constructional  definition 
s 


dx-,  .  . .  dXj^ 


(4.23)   n^(^iPi  ...  ^3^3)  =f     ^Wj    ^^1 


dx,   .  . .  dXj,    exp  -  [   ] 
^1       ^M-l 


t  I         It 

p(x^P^)  ...  p(XgPg) 


where 


N 

(4.24)  pi^^\)=Xl^i^±^   -  ^k)  • 

Since  the  external  potential  (4.l6)  is  expressible  as 

M 

(4.25)  U  =  211  e  j  dXj^  P^^k^k^  ^(^A^  ' 
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the  truth  of  (4.21)  Is  easily  verified.   Formulas  like 
(4.13),  (4.1?)  -  (4.20)  bear  a  resemblance  to  the  corresponding 
expressions  for  classical  mixtures,  with  the  temperature 
index  taking  the  place  of  species  parameter.   The  inter- 
actions from  the  mixture  point  of  view  are  somewhat  peculiar, 
see  (4.12),  in  that  one  of  them  involves  particles  of  the 
same  species  only,  while  the  other  is  nearest  neighbor  in 
species  index,  see  (4.11).   The  resemblance  is  exploited 
in  the  final  two  sections.   See  diagram. 

^K-,  ^K         ^i/^K  +  l 


AX. 


/I 
/  \ 

/    V2 


/ 


I 


I 
I 

I  ' 
I  I 
\  > 

\  V 
\  \ 
\  \ 


^2 


^I'^N 


'"       X, 


^OX; 


N 


^N 


The  interactions  in  the  exponent  of  the  path  integral 
(4.12). 

While  we  have  so  far  taken  expectations  of  coordinate 
operators  only  (6-functions ) ,  the  path  integral  representation 
easily  accomodates  momentum  expectations.   These  will  involve 


21 


full  self  distributions  at  inf initesimally  separated  temper- 
ature arguments.   Thus  for  example. 


(4.26) 


N   .   .  N 

im  •     ^    m 


i=l  1i   ^a      he    f^       ip     1r 

^k  ^k     "^k-l 


=  >  Z^m^      j  dx^  ..  dx^J   dx   ...  dx 

^1       Pm-1 

N 

im 


exp  -  [   ]  '  XZ   ^T  (^n    -  ^n       ) 


^k    Pk-1 


(4.27)       =  11m   i^  /  dxdy  n„(xp,yp,  t)(x  -  y) 

e_>  0  ne  ^       ^   k.  k-j. 


which  is  of  course  zero.   One  must  be  more  careful  with  the 
kinetic  energy  expectation.   Feynman  shows  [8]  that 

N 
(4.28)    -  ^~-r^  (^i    -  ^i  Hx.  -X    ) 


has  a  finite  expectation.   This  works  out  to 


(4.29)   <  K.  E.  >  =  -  lim   — 5-^  /  dxdydz(x  -  y)(y  -  z) 

e->0  2h  e"^  J 


'  ^3^^^k+i  ypk  ^Pk-i)  • 
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5-   Classical  Statistics:  Bootstrap  Operation 

From  here  on  the  Integral  representation  will  be  used 
exclusively.  Formulas  (4.13),  (^.17),  (4.21)  -  (4.23)  are 
basic . 

In  this  section  we  derive  relations  between  lower  and 
higher  distributions.   The  general  idea  is  that  a  distribu- 
tion defined  for  a  system  in  which  one  or  more  particles  are 
"fixed"  is  convertible  by  integrations  into  a  higher  order 
distribution.   To  "fix"  a  particle  means  to  omit  integrations 
over  it  at  every  point  of  its  path.   Most  of  these  integra- 
tions are  restored  by  the  conversion  process  mentioned. 
This  is  a  formal  procedure.   There  are  no  recoil  effects, 
the  kinetic  energy  acting  like  a  constant  throughout. 

By  fixing  a  number  of  particles  in  the  partition 
function  all  higher  distributions  may  be  obtained.   Such 
a  relation  is  the  quantum  analog  of  (2.15). 

Let  us  write  the  partition  function  for  a  system 
in  which  j  particles  are  fixed.  These  provide  an  external 
potential  for  the  remaining  N  -  j  particles,  thereby  intro- 
ducing a  non-uniformity  intoq. 
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N-j 


(5.1)  ^[u^]  =ZlTNr37T  J  d^j+i  •••  d^NJ  --J+1.  •••  -N 


N=j 


(ix  .  , ,    ...    dx,T   /    dx  . ,  T      ...    dx. 


■Pn 


P- 


M-1 


M 


exp   -    e 


^     m      -2      ,    T^       m 


S_^ 


X  . 


k=l    ^  1=1   2h^      -^a     l=j+l   2^2      ^P 


k 


1      J  1      ^ 

+  ?  IZ  Hx.         -X         )+2XlMx.       -X         ) 

1  Pic  ^k  J+/  ^k        ^k- 


M        N 


exp   - 


e^E 


.<c 


k=l    i-j+1 


^k 


(5.2)  U^(x        ,P^)    =|ZMXi        -  x^^    ,p^)    , 


■Pk  ^Pk 


[cf.    (2.13)].      For   convenience   the   constant    first   and  third 
terms   In  the   exponential   have   been   Included.      Now  perform 


(5. 5) 


•J   = 


J 


dx-,    .  .    dx  .    /    dx- 


■P- 


/   dx  . 


'M-1 


I  t 

5(x-,      -  X,  )5(x-,       -  x^) 


5(x 


^n 


^A+l 


5(x.      -   X    )    , 


each  6  eliminating  an  Integration.   The  canonical  part 
transforms,  by  (4.17): 
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(5.4)   I^   -(N-j)r-   ^  "^s    (^1^1  ••  ^n^n^  ^n+l^n+l  "    ''    '•- 


t  t   ,7 
s^s   NT 


N 


N         ^ 
Multiply  through  by  z  and  sum  ^ 


fe 


(5.5)   Z''lJ^[U|,)  =  n^J>(x;p;  ..  ^X;  ^n+lCl 


I  t 
.  .  .  ;  ...  ^gPg ) ^  ' 

This  is  the  analog  of  (2.15).   The  final  distribution  which 
we  get  on  the  right  of  (5.5)  is,  of  course,  for  a  uniform 
system. 

Some  special  cases  of  (5-5): 

(a)  J  =  s,  which  means  one  integration  is  stopped  per 
particle, 

(5.6)  z"^  /  dx-,  ..  dx  .  /  dx,   ...  /  dx  .    6(x,   -x,  )5(x„   -x„) 

Pi       pm-i   Pi      P2 

J 

a  completely  mutual  distribution. 

(b)  J  =  1,  one  particle  fixed, 

(5.7)  zj    dx^y  dx    ...  dx     5(x    -x{)...5(x    -x;)j.[U^] 

^1       ^M-l    ^1  ^s 

t  t  t  I       II 


=  ^^(^1^1X2^2  ••♦  ^s^sM' 
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a   full   self   distribution. 

(c)      with   s    =  2   in   (5.6)    and    (5.7),    we   have   formulas 
for  the   mutual   and   self  pair  distribution  respectively. 

Suppose  we   had  begun  with,    in  place   of    (5.I), 


(5.8)      n|^^(x"p]^'     . 


II      II  II        ti 

• •    ^n  ^n    ^    ^n+l^n+1        ' ' 


\h    ^''^y^^ 


-T 


N-j+k 


(N-j-k)] 


d^j+1    ..•    dxj^j^    ...    dxj^ 


dx  . ,  -,    ...    dx^T 


exp 


M-1 


^    PJ-     m      .2 


Pi. 


N 


+    >  TT    X  . 


x^     +  iZZ  Mx.    -  X.    ) 


k 


±M 


Pk        ^Pk 


1     ^ 
+  p  ^  Mx.       -   X.      ) 


i^i 


k 


P 


M        N 
exp-e  >       y~ 


k 


Ui(x        ,p    ) 
k=r  T^+1      'P      ^p^     ^ 


5(x      -^        -  x^    ) 


"P 


n+1 


ti 


Pt 


where   the   arguments   l,...,t   are   consecutively  distributed 
among  k  groups    set   off  by  j_  (particles    j   +1,...    j   +  k),    and 
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with  particles  1  ..  j  fixed,  supplying  Ui,  as  before. 
Then  If  (5.3)  Is  applied  to  (5.8), 

/  ^  i     1       f  k)  ,     "      "  II      II  ti  II 

(5.9)       zJlJ    n^^^x,    P^    ...    x^   P^    ;    x^^,    P^^^    ...;     ...; 


II      II 

X 


t    h    lU^^pfU^^ 


"^sit      (^1^1    •••    ^n^n'    ^n+l^n+l    •"'    '"'    "'    ^s   ^s    ' 


II      II  II      n  II         II  ti      11 

^1   Pi    •••    ^n  Pn   ^    ^n+l^n+r--    ---;•••    ^^   ^^    )^     , 


which  Is  the  generalization  of  (2.14). 

Example:   j  =  l,  s=l,  k  =  l  and  t  ==  1, 

(5.10)   zy^dx^y  dx^^  ...  dx^^   5(x,^^-  x;)n,{x"p;'  |U^)^[U^] 


=  ngCx^P^;  x^  p^  )p 


the  mutual  pair  distribution. 

In  addition  to  (5.5)  and  (5. 10)  separate  formulas  can 


be  derived  for  the  n  . 

s 
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To  see  this  write  the  canonical  partition  function  for  a 
system  of  N  particles  In  an  external  potential: 


(5.11)   Z^[U]  =  /  dx^  ...  dXj^  /  dx-j_   ...  dXj^ 

Pi       Pm-1 


^M   r  .1   _   r.  N 

exp  -  e  ' 


Pk  J^i  ^    Pk 


J    N  N  - 

1  j+i     'Pk    \      1       Pk     • 


(5.12)  =  f  dx,  ...  dx  .  r  dx-,   ...  dx.    Z^'-^'CU  +  U.] 


fa-  J-  1   2ii2   ip,^    2  ^  ip 


x^   ) 
k       -         ^k    ^k 


+ 


d^k  C(Vk)U(x^p^)]  , 


(5.15)  C(^A)  =  SI  5(^1^  -  ^k)  ' 

1     Pi 


k 


with  Ui  as  In  (5-2). 

Differentiating  both  sides  s  times  with  respect  to  U, 
using  (4.21), 


oft 


(5.14)   n^(x^P-^ 


x^P^|U)zN[U]   s. 


dx-j^  ...  dx  . 


dx^   ...  dx  .     exp  [   ]^^ 


!   I 


t   I 


<^N-j 


I    I 


t  t 


I^C(x^p^)  ...  ct^tPt^'^sIt^^t+A+r--  \.Pslu  +  u^) 


.  7.N-J 


'•      9-'  (s-t )  It  i 


Taking  s  =  j  =  1  for  example. 


(5.15)   nN(x|p{|U)zN[U] 


dx-j^  /  dx,   . . .  dx,     exp 
^1       ^M-1 


-  e 


'^   ~  m   -2 


p, 


+U(x    ,p^) 
^k    - 


Ux'p')zN-l[u  +  Ui] 


'1^1 


>■ 


^5^-i(xXlu  +  u^)zN-i[u  +  u^]]  , 


this,  by  (4.23),  (4.24)  and  (5-13), 


dx,  /  dx. 


dx. 


Pi 


■P, 


exp  -  [   ] 


U 


M-1 


N6(x    -  x;)zN-l[u  +  U^] 


which  goes  over  into  grand  canonical  form 
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I  t 


(5.16)   n-^{xjp-[|U)^[U]  -  z  /  dx-j^  /  dx^^  ...  dx^^   exp  [   ]y 


■Pi 


"P, 


M-1 


S(^l^,-^i;^fU+U^]  , 


■P 


which,  except  for  U,  Is  (5-5)  with  j  =  s  =  1. 
With  J  =  1,  s  =  2 


(5.17)  nl{x[^y^^'^\\J)Z^[\J] 


j^^lj 


dx-,   . . .  dx,     exp  -  [   ] 
^1       ^M-l 


U 


[ 


C(x|p|)C(x2P2)2^"'^tU  +  U^] 


'  '  ,aN-1, 


I  I 


N-1, 


+  C(Xl6l)n5;--'(x2P2lU  +  U^)Z^^-"[U  +  U^]  '2 


aN-1 


I   I   >   ! 


N-1, 


+  ;\^-^(x^P^X2P2lU  +  U^)Z^^-^[U  +  U^] 


dx. 


exp  -  [   ]. 


^    ^1      Pm-1 


N6(x,   -  x')5(Xt   -  x')Z^  ^[U  +  Ui] 


■P^ 


^p; 


' xaN-1, 


t  I 


.N-1 


+  N6(x-L   -  X^)n5|--'(X2P2|U  +  U^)Z    [U  +  U^] 
Pi 
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till 


(5.18)   n^(x-Lp^X2P2lU)^[U] 


z  /  dx-j^  /  dx-j^   . . .  dx-,     exp  -  [   ]y 
^1       ^M-1 


[6(x    -x;)5(x    -x;)^[U+U^] 
Pi  ^2 


on  the  grand  ensemble. 

To  summarize,  we  have  in  this  section  found  relations 
between  lower  order  distribution  with  external  potentials 
(due  to  fixed  particles)  and  higher  distributions  for  uniform 
systems.   The  main  results  are  embodied  in  (5.5),  (5-9) 
and  (5.14).   These  relations  will  now  be  applied  to  the 
problem  of  obtaining  Integral  equations  and  bounds  for 
distributions . 


6.   Bounds  on  Distributions 

We  shall  require  in  this  and  later  sections  a  functional 
Taylor  expansion.   Consider  a  function  of  space  and  tempera- 
ture arguments  a)(xp)  which  changes  from  initial  to  final 
value  due  to  a  change  in  a  parameter  a,  a  varying  from  0  to  1; 


(6.1)   a3(xp,a)  =  co^(xp)  +  a[(D(xp)  -a)Q(xp)] 
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Given  a  functional  of  (Xi(xp),  ^[a)(xp)],  we  wish  to  expand  it 
in  the  deviation  of  o)  from  its  initial  value.   Since  f   can 
be  regarded  as  a  function  of  a  through  (6.1),  the  Taylor 
expansion  for  it,  including  remainder  is  just: 


(6.2)   ^(a) 
a=l 


^ —  J 


J-0 


ha' 


a  =  0 


+  /  '  i!!^  (l-f'  da 


Since  by  the  chain  rule 


(6.3) 


h^r  ^co(xP,a)     6    dxdp  , 
Sa   /    5a     5(D(xP,a) 


(6.4)      =  /  Aa)(xP)  s^.(.;^^^)  dxdp  ,   (Am  =  CO  -  cd^)  , 


the  functional  expansion  is 


(6.5)   ^[cd(xP)]  -^[co  (xP)]  +  XI  TT  /•••[  A^(^lPl)---  Aa>(x  P  ) 


J-1 


6a3(x-]_P-|_)  .  .  .bcDix  -P  •) 


CD    =    CD 


dx-j^dp^...    dx.dp. 


o  ^' 


Acd(x^P^)...    Aa3(Xg^^pg_^-j^) 


6^+^^[cD(a)] 


5CD 


,    Q,     ^10^^"^  dx.dp,  ...dx_-,dp_.da 


For  example,    if  cd(xP)    =  U(xP),    f[(ii]    =^[U],    and  cd^  -  0,    we 
have 
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{6.6)     -xiu]   -;^[0]  +  ZI -^T^  J   ^i^Pi  ••  ^^j'^Pj^f^^ 


p 


+  J        J    dx^dp^  ..  d^3^idp^^^da  i-i-g^fL-  (_l) 


.[0]n^^^(x^P^  ..  ^3+1^3+1, a)U{x^p^)..  UCx^^.P^^-^) 


In  which  (4.21)  has  been  used. 

Of  particular  Interest  is  the  case  of  purely  repulsive 
potential  U.   Then,  since  the  n.  are  positive,  the  Jth  term 

J 

in  the  series  has  the  sign  (-1)"^,  while  the  remainder  has 

s+l 
the  sign  (-1)    .   The  series  alternates  and  successive 

upper  and  lower  bounds  may  be  iterated  out: 


(6.7)  ^[U]  <  ^[0] 

>^[0]  -J  dx^dP^n(x3^P^)p.[0]U(x-LP^) 

<^[0]  -  /  dx-^dp^n(x^p-^)^[0]U(x^P^) 

If  we  specialize  the  U  to  the  pair  potential  arising  from  t 
fixed  particles  at  y,  . .  y,  , 

t 

(6.8)    u.  (x  p, )  =  XI  M^i  -  y.  .Pi)  > 
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bounds  on  uniform  distributions  may  be  obtained.   For 
example,  if  t  =  1,  by  (5.16), 


AO, 


(6.9)   ni(yiPi)  in^iy^^^) 

>n°(y3_P^)  -  /  di^dJ^dP^H^Cx-^P-^) 


^o, 


A,  o 


or  with  t  =  1  and  t  =  2,  by  (5-5)  and  (5-l6), 


AQ, 


(6.10)  rigCy^P^ygPg)  1  '^2^yih^2^2^ 


>  ng  - 


J   d^dx^dF-Ln°(x^P-L) 


O/'^-r 


(t)(x-^  -  x)  n^CxP^y-^P^ygPg) 


, -^o, 


+  n°(xP-Ly-LP-j^)K^(y2P2) 


^o 


+  n^(y-LP^)n^(xp^y2P2) 


(The  Interaction  between  fixed  particles  [t  =  2]  has  been 
neglected. )   This  is  the  beginning  of  a  Mayer-Montroll  type 
sequence.   Once  the  n  are  known  the  individual  distributions 
may  be  found  by  the  equations  of  Section  5. 
If  {6.6)    is  differentiated  t  times. 
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I   !  t   I 


(6.11)   n^(x-Lpi  ..  \^t\^)^[^]    = 


•  ^t+T^^l^l  ••  ^t^t'  ^1^1  ••  ^j^i^ 


•  U(x3_P-l)  ..  U(x  p  ) 


J  J 
:^dp^  ..  dx^^,dp^^,da(iz£).' 


It      It 


>[o]^+s+i(^iPi  ••  ^t^t'  ^A  ••  ^s+A+i'^) 


u(^iPi)  ••  ^(x^^iP^^-l)  , 


and,  for  example,  choosing  t  =  1, 


(6.12)  n^(x^P^|U)^[U]  <  n^(x-^P^)a^[0] 

-  f  dx^dp^[0]n2(xjp|x^P-L)U(x^p-L) 
Specializing  U  by  fixing  one  particle, 

(6.13)  U,(x^P^)  =  H^i   -  Vi      >^i)> 

Pi 
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and  using  (5.I8),  [also  (6.9)], 


I  t 


'-O 


(6.14)   n2(x-LP^y-j_P-|_)  <  n^ 


>  n°  -  /  dxdx-,dP-,  <t)(x,  -  x) 


^  O  /    o      \     O/^TT    n  >  'vO 


n^(x^P^)n^(xp^y-LP^)n^(x^P^) 


+  n°(xp^x|p|y^p^)n°(x-^p^) 
+  n2(x-j^P-|^x^P^)n2(xP^y-]^P-|_) 


This  with  (6.9)  is  the  beginning  of  a  Klrkwood-Salsburg  type 
sequence. 


7.   Classical  Statistics;  Debye-Huckel  Sequence; 

Direct  Correlation, 

In  Classical  Statistical  Mechanics  It  Is  possible  to 
derive  a  general,  self-consistent  sequence  of  Integral  equations 
for  the  modified  Ursell  functions  F^  [7].   The  simplest  of 
these  is  the  linearized  Debye-Huckel  equation  (2.17). 
Unfortunately  this  procedure  cannot  be  generalized  to  quantxam 
mechanics  because  the  quantum  F^,  (4.22),  do  not  satisfy 
bootstrap  operations  of  the  kind  developed  in  Section  5. 
We  therefore  confine  ourselves  to  the  pair  distribution. 
We  seek  a  consistent  approximation  for  the  self  and  mutual 
distributions. 
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Expand  log^[Ui]   v.  Ui  where  r[U,]  given  by  (5.I) 
J  =  1  and 


(7.1) 


U^(^i   ,Pl,)  ==  H^.      -  X        ,^^)    , 


■P 


k 


\       ^\ 


particle  "1"  fixed. 

The  functional  Taylor  expansion  Is 


(7.2)   log^[U^] 


r  P 5  log  ^[Ui] 

=  log  a  [0]  +  /    dPdx  1_  'P 

'^      ^0  6U^(xP) 


Ui  =0 


^^ 


U^(xP) 


p 


+ 


dPdxdpdx 


^^g^^Uj] 


0 


6U,(xP)6U,(xP) 


>^Kxp)ou^ 


Ui  =  0 


'^ 


(7.3) 


U^(xP)U^{xP)  +  ... 


P  -,^ 


=  log^[0]  -  y  dpdx  n-j^(xp)<t)(x  -  x^  ,  p) 


P 

+  I  /    dpdxdpdxfn^CxP;  xp)  +  n^Cx^xp") 

-  n^(xP)n-^(xp)] 
.  (t)(x  -  X,  ,  p)(t)(x  -  X,  ,p)  , 

p  p 

In  which  use  has  been  made  of  (4.22)  and  (3.I8). 
Differentiating  (7.3): 
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+  n^lx"  p"xp)  -  n-^lz"  p"  )n^(xp)) 


.  ([)(x  -  x-j^  ,P)  . 

F 

Linearizing  (7-3)  and  dropping  an  inessential  term  linear  in 

(7.5)  ^[Ui]  =^[0]  +1/    dpdxdpdi?^[0] 

(n^tnp;    xP)    +  n^CxPxp) 

-  n^)(i){3r  -  X-,    ,F)M^   -  ^1    .P)    • 
—  & 


Applying    (5-7) :    (J    =  2,4) 

^P ^  ^ 

(7.6)      ngtx'P'x-'P")    =  n°(x'P'x"  P"  )    +  |j    dpdxdpdxdxdx 

0 
(nglxP;    xp)    +  ngCxPxP)    -  n^) 


<t)(x    -  x)(|)(x    -  x) 

•    nj(x'p'x"  p«'Srpx"^)    , 


the  n?   being  a  completely  self  k   particle  distribution, 
"o"  denotes  free  particles. 

Applying  (5-10)  and  (5-7)  (with  j  -  1,2)  to  (7-4), 
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after  one  has  multiplied  through  by  z^[Ui ], 

P  _  _ 
(7-7)   n^Cx'p';  x"  P"  )  =  n^  -  /  dpdxdx(n2{x"  p"  ;  xp) 

0 

+  n2{x"  P"x3)  -  n^)<j)(x  -  x) 

n^Cx'P'xF)  . 

Equations  (7-6),  (7-7)  are  consistent  simultaneous  integral 
equations  for  np(  )  and  n^{;),    valid  for  sufficiently  weak 
potential. 

A  cruder  approximation,  which  has  the  virtue  of  being 
exactly  solvable,  is  to  replace  n^{  )  by  its  free  value  in 
(7.7): 

(7.8)  n(x'p';  x"  P"  )  =   -    f  dpdxdx(n(x"  P"  ;  x^) 

+  n°(x"  P"xP))(t)(x  -  x) 
•  n°(x'p'xp)  , 
n(; )  =  n^C;)  -  n^. 

The  solution  of  (7.8)  proceeds  as  follows. 
The  Fourier  transform  of  (7.8), 

(7.9)  n(k,P'P"  )  =   -    f     dF(n(k,P"p}  +  n°(k,P"p))<t)(k)n°{k,pp') 
has  for  kernel  of  the  homogeneous  equation 


^ 


(7.10)   n°(k,pp')  =  n  exp  -  ^^  (|p'  -  pi  -  l^''^'^  )  , 


see  (3.25). 

This  kernel  has  the  property  that 

(7.11a)   n°(k,P  -  a)  =  n°(k,a)  , 

a  =  IP'  -  Fl 
and  we  define,  for  temperatures  >  p, 


(7.11b)   n°(k,p  +  a)  -  n°(k,a) 


It  then  follows  easily  that 

(7.12)   n(k,P'P"  )  =  ;^^-n_  ^^(pn  )^^(p,)  ^ 

where 


n     n 


(n  =  +  int,  0. ) 


O/,  „x  27rlna 


(7.15)     \{k)    =    1        da  n^(k,a)  exp  £1^ 

(7.14)      ^^O')  =  p-1/2  g^p  _  27^  ^ 
and 

The  A^  and  ^^   being  eigenvalues  and  elgenfunctions  of  the 
kernel.   Once  .the  pair  distribution  is  known  one  can  work 
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backwards   to  find    o^and  hence   the  thermodynamics.      This   1: 
done   by  recognizing  that   if 


^  — >  7<t>      , 


0    <  7    <   1 


(7.16) 


d   log  ,-^7)    ^ 


-  f     <  7^   >y    , 


(7.17)  log  ^(1)    =  log  3,°   -   p   /        ^     <7<t)> 


7 


Since,    setting   P'    =  P"     in    (7.12) 


(7.18)  «t)>  -  I  (27r)-^  V  /   dk  n      fk)    <t)(k)    , 

^  (mutual) 


we  arrive  at,  by  (7-12),  (7-17)  and  (7-18), 


o  .  1 


-5 


(7.19)   log  ^  =  log  ^^  +  i  (27r)~''  V  /  dkj 


Mk)A 


n   — 


n 


-  log(l  +  H^)\) 


This  is  the  result  obtained  by  Montroll  and  Ward  in  their 
theory  of  the  election  gas  [10]. 

At  high  temperatures,  p  — >  0, 


n 


^iP^o  ' 


(7.20) 


n(k)  =  - 


n^pMk) 
H-nP(|)(k) 
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which  Is  also  the  solution  of  (2.17),  and  gives  rise  to 
the  Debye-Huckel  equation  of  state  for  a  screened  Coulomb 
potential  (neutralizing  background)  [2],  [10]. 

Equation  (7.8)  Is  the  quantum  generalization  of 

(2.17). 

Returning  to  classical  statistical  mechanics  momentarily 

we  recall  that 


(7.21) 


5n(x 


-5pU 


U) 


yj-  =  n(x|U)&(x  -  y) 

+  n^(x,y|U)  -  n(x|U)n(y|U)  , 


is  composed  of  a  direct  effect  at  y  due  to  U(y)  plus  a 
longer  range  reaction.   Let  us  consider  the  inverse  derivative 
(which  will  have  meaning  on  a  grand  ensemble  where  arbitrary 
density  variations  can  be  made)  and  split  it  into  singular 
and  non-local  parts: 

^^■^^^  6nUtu)  -  n(x|U)    -c(x,y|U)  . 

c(x,y|U)  is  called  the  direct  correlation  function  [1],  the 
name  to  become  clearer  presently.   The  derivatives  in  (7.21) 
and  (7.22)   are  matrix  inverses  of  one  another  and  therefore 


(7.23) 


5n(x 


-5PU 


U)   -5pU(y) 

rn  6n(ztu)  ^y 


5(x  -  z)  . 


Substituting  (7-21)  and  (7.22)  into  (7.25) 
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(7.24)    g(x,y)  -  1  =  c(x,y)  +  /  c (x, z)n( z ) (g( z ,y )  -  l)dz  , 


where 

np(x,y) 


(the  presence  of  U  is  understood). 

Equation  (7-24)  when  Iterated  becomes 


(7.26)  g(x,y)  -  1  =  c(x,y)  +  /  c(x,z)n(z)c(z,y)dz 

+  /  c(x,z)n(z)c(z,oD)n(cD)c(a),y)dzda) 

I.e.  g  -  1  Is  represented  by  a  "direct  correlation"  from 
X  to  y  plus  a  siom  over  all  the  ways  direct  correlations  can 
be  transmitted  via  Intermediate  particles  from  x  to  y  at 
density  n. 

One  expects  that  c  will  have  a  shorter  range  than  g  -  1, 
of  the  order  of  the  Interpartlcle  potential.   (The  non- 
Interactlng  g  Is  unity).   In  fact,  for  hard  cores  c  has 
exactly  the  range  of  the  potential  [9]. 

Consider  again  the  Debye-Huckel  equation  (2.I7).   Here 

(7.27)  (g  -  l)(k)  =  -  ,itW,   .   [7. 20] 


In  general,  from  (7.2^),  U  =  0, 


(7.28)         (g.i,(k)  =^^£l|^. 


itj 


so  that  In  this  case 

(7.28')  c(k)  -  -  p<|)(k)  , 

the  potential  Itself.   Knowledge  of  c  supplies  the  thermo- 
dynamics, for  as  may  be  shown  [11]  (U  =  0), 


(7.29)  p  ^  =  1  -  n  j  c(x)dx  , 

p  =  grand  canonical  pressixre,  =  log  ^(VP)~ 


In  an  attempt  to  generalize  these  ideas  to  quantum  systems 
we  have  first,  in  matrix  notation. 


^^•^°)     -:^  =  ^2^      ^  "^  ^2^'^  "  "^  ' 


((4.22)  with  s  -  2) 


The  inverse  derivative 


(7. 51)  -  |§  =   n-^   )  -  c  , 


with  n~  the  matrix  inverse  of  the  self  distribution,  which 

5  f  X  '  -X  "  ) 
at  equal  temperatures  is  ^^ ,  and 


(7  32)  _^  .  z5U  -  1 

^l'^"^'  -5U     Bn  -  -^  • 
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Substituting  {7  .J>0)    and  (7.31)  into  (7.32)  we  have  for 
uniform  systems  (U  ~  0) , 


(7.53)  n  =  n  c  ng  +  ng  c  ng  , 

n  =  n^  ( ;  )  -  n 
It  is  convenient  to  replace  c  by  n~   c '  nl   ,  leading  to; 


(7.3^)  n  -  c  +  n  n^^  c  , 


dropping  the  prime. 

This  iterates  in  the  fashion  of  (7.26) 


(7- 55)  n=c+cnp     c+cn^     cn~     c+ 


The  thermodynamic  relation  can  be  shown  to  be 

A  /  P      c  ( X ,  a ) 

(7.56)       P  ^  =  1  -  I  /    dxda  2—   • 

^^    J  Q  n 

As  an  illustration,  consider  (7.8).   There 


(7.37)  n°(k,P'P")  ^XZK   K^^''  ^K^^'^    ' 

n 
so  that 

-1  -,      ^ 

(7.38)  n°   (k,p'P")=IZ  A"^  ^^(P"  )<t)n(P')  , 

n 
and  ry 

(7.39)  n(k,P'P"  )  =YZt^  Ki^''  ^K^^'^    • 

n     n 
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The  Fourier  transformed  equation  (7.3^)  is 

(7.40)   n(k,P'p"  )  =  c(k,P'P"  )  +  /    n(k,P'F)n:-^(k,FP) 

^0 

•  c(k,pp"  )dpdp  . 


Making  the  substitutions  (7-38)  and  (7.39)  one  Identifies 
(7.41)    c(k,P'P")  =  ^^  <t,*(pn  )<|,^(p,)  . 
The  classical  limit  of  this  Is  of  course  (7.28'). 


8.   Quantum  Statistics:  Toron  Decomposition; 

Bootstrap  Operation 

In  this  section  relations  like  those  of  section  five 
will  be  developed  for  the  Elnsteln-Bose  and  Ferml-Dlrac 
cases.   We  confine  ourselves  to  one-  and  two-particle  distri- 
butions only.   We  wish  to  relate,  for  example,  the  one- 
particle  distribution  in  external  potential  due  to  fixed 
particles  to  the  pair  distribution  for  a  uniform  system. 

Before  this  can  be  done  it  is  necessary  to  rewrite^  and 
the  distributions  in  forms  which  display  the  permutatlonal 
sum  in  a  way  convenient  to  our  purposes.   These  new  forms 
also  permit  easy  calculation  of  free  particle  functions. 
Recall  that 
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(8.1) 


,N 


^ 


+  ^P 


dx-|  .  .  dx 


N 


1q      N 


P 


M-1 


M 


r  N 


exp  -  Y. 


m 


far     -  f-r  2t2 


(Normalization  factors  are  omitted  for  convenience  in 
writing.)   Consider 


(8.2) 


+  ^P 


z,-2Z(t) 


dx,  .  .  dx-j  /   dx-,   . .  dx^     •  exp  -  [  ] 

^P   Pi     Pm-1 


From  these  permutations,  we  extract  a  "closed  loop"  of  t 
particles  interacting  with  the  remainder,  t  ^  N 


.oop  =  ® 


closed  loop  =  (j)   dx,    . .  dx, 

P^       P 


M-1 


means   x,   =  x,   ,   x„  =  x-,   ,  .  •  • 

^o    ^P     '^o  -^p^ 


X,   =  X,  -,  ,   i.e.  particle  one  returns  to  its 
^o    "^"-^p 

starting  place  after  "time"  tp.   In  the  language  of  Montroll 
and  Ward  [10]  this  is  a  t  particle  toron.   There  are 
N(N-l)  ...  (N-(t-l))  •  i-  ways  to  construct  such  a  loop,  the 
i-  factor  eliminating  redundancies.   This  much  of  Z^^  now 


appears  as 


^7 


(8.5)  2Z  (-)^  f  <^^i  •••  d^N  (-^^^"^  M^i   •••  dx^ 

P      ^  ^    Pi      Pm-1 

^'    """  -  '   dx,  ^T   .  .  .  dx,,     exp  -  [  ]  , 


(N-t)!   N   J   -t+1   •••  -N 

the  permutation  sum  operating  now  only  on  particles  t  +  1,...,N, 

The  end  points  are  as  usual  x^  =  x, ,  .  . . ,  x^  =  x^,  .  .  .  ;Xj^  =  x^^. 

o  o  o 

The  full  Z,,  Includes  a  sum  over  all  sizes  of  torons 
N 

t  =  1, . . .  N.   So  that 

CO  ^N   N   /+Nt+1  ^ 


dx^T  <P  dx-,    . . .  dx. 


J    dx^  ...  d^^j   ^^^   ...  ^^^ 


^1        ^M-l 


^t+1^  •  •  •  ^^N^    exp  -  [  ]  . 


Reversing  the  order  of  the  N  and  t  sums  with  N-t  — >  s, 
we  have  finally 

OO    00        r7S+t        ,  ,T   


;  dy,  ...  dy.  ^  dy,    ...  dy. 
J    1       ^^   ^Pt       ^P, 


^1        ^M-1 


j  dx^  .  .  .  dXg  j_  ^^^^  •  •  •  -s 


dx-,   .  .  .  dx 
P  \  'Pm-1 
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exp  -  e 


M   r- 


T-^  m   -2   ^  v^  m   -2    ,1  ^^— 'i, 
2ll-T^y±     +  X_  -r^  ^1   +  n-  y^  Mx   -X.   ) 
^r   2t^   ^P,^  ^T-   di2   ip     .f  ^-    1     J 


k 


"\   Pk 


t    s 


+  1^  '  Myi  -  y.  )  +X:±:  My, 

1  3,.     -^a,.        11        ^ 


'Pk     ^Pk 


Pk     '\- 


The  third  term  in  the  potential  is  the  toron  Interaction 
with  the  rest  of  the  system.   The  second  term  is  the  toron 
self  interaction. 

In  a  similar  fashion  this  toron  decomposition  may  be 
accomplished  for  n,,  n^{    )  and  n^C;)  or  alternatively 
these  expressions  may  tie  derived  from  (8.5)  by  differentia- 
tion with  respect  to  U. 

We  record  the  results: 


-1 


Z^^^  -^-  ,+^p+t+l 


(8.6)   n^{x'P')  =>-^T-2Z  V5I(^) 

t=l  s=0       p 


J    ay^    ...  dy^  j^  dy    ...  dy^ 


P. 


'P 


M-1 


dx,  . . .  dx   /   dx^   ...  dx     exp-[  ]  , 
^     "^1       ^M-1 


-1 


5(yi  -  xM  , 
■Lp« 


,s+t 


(8.7)   n2(x'P'x"  P"  )  =jr     ^Z  2Z  ^^sT-  2Z  (-) 


+  .p+t+l 


t=l  s=0  "•    p 


exp 


-  [  ]  &(y-,  -  xM&Cyn   -X"  ) 

-Lpi  -LptI 
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and 


(8.8)   n„(x'P';x"P"  )  =  J?J^  J~  y~ 
^  ^      t=l  s=l 


F7~  2 —  ^-' 


exp  -  [  ]'6(y,   -  x')5(x,    -x "  ) 


^>'^^^ZI(t)P^^^^ 


X^  i^  "•  T 


r  ...  f  exp  -  [  ]  •  XI  5(yn   -  x')5(y.    -  x" 

Formulas  (8.6)  -  (8.8)  give  a  relatively  easy  way  of  calcula- 
ting free  distributions. 


(8.9)   n°(x'P')  =;^  (1)^''^  z*  /  dy- 


^^t 


i  dy^   . . .  dy^    *  5(y-j_   -  x') 

^    Pi     Pm-1     P' 


1^   t   ^   ^ 
exp  -  e  ^__  ^_  -j-^  y . 

k=l  1=1  2Ti    Pj^ 


but  by  (4.10)  this  is 


(8.10)   n°(x'p')  =-f-   (j:)^+^  z*  jK°(y^Py^O)K°(y^Py^_^0) 

...K°(y23x'P')K°(x'P'y-L0) 


dyn  ...  dy.  , 
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(8.11) 


t=i 


t 


3/2 


by  (3.7)  or  (3. 8), 


Since 


A  .  (-4L_)V2 


_  z  a  log  -^- 
1  ~  V  ^ — 


(8.12) 


log  3.°  =  vy-   (+)t+l  _z   ^ 

"       n  +--5/2 


Similarly 

(8.13)   n°(x'P'x"p")  =  ^T-  ^+^t+l  „t  / 


^(i)^^^^V^^i 


dy. 


O'  >  P"  ) 


K°(y;^Py^O)K°{y^Py^_^0)  ... 


K°(yaPx'P')K°(x«P'x"  p"  )K°(x"  P"y^O)  , 


(8.14)   n°(k,a)  =  y-  (+) 


^     -aa  +aa 
372   ^® 


a  = 


+-2,  2 

^^  -  =  IP-  -  P"  I  , 


and 


(8.15)   n°(x'p';x"p"  )  -  n^  =  J^ 


t-2 


(+)t+l  ^t 


1=2 


dy-^  . 


•  dy. 
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rO, 


0>  >  P"  )     .  K^Cy^py^O)  ...  K"(y^^^Px'P') 


•  K°{x'P'yiO)  ...  K^Cy^Px-'P"  ) 


•  K^(x"P"yiO)  , 


(8.16)   K^  (k,a) 


,  -    t    t-1    -aO^+a)+a  lP|i^ 
t=2         t-'/'^    i-1 


As  a  check  on  (8.I6)  transform  back  into  x    space  at  a  =  0; 


00     t-1 


(8.17)   n°  (x,o)  =  T~y~  (+) 


? 


t+1 


zH2 


2  I3T 


t3/2(,.  |,3/2 


^    t 


Redefining  Indices  t    -   £   — >  m. 


00      00         „ni+^  vr,_l_/)j_1      O 

(8.18)  E°(:c,0)  =21:11-^-172  (-)     ^ 


exp  [-  x^TTA^/^Ci"-^  +  m"^)] 


+  .2 


r  00   „m     ^ 


-  m=l  m'^' 


the  familiar  London-Placzek  formula. 

One  other  free  particle  function  will  be  included 
here  for  future  reference:   the  free  modified  Ursell 

Aq 

functions,  F  , 
'   s' 
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(8.19)  F°(x-,P-,  .. .  X  p  ) 


oo 


t+1  t 


y~  (i)     z   /  dy-j_  .  .  .  dy^  c6  dy 


^1 


dy, 


'P. 


M-l 


M    t 

exp  -  s  n:  H  ^  y? 

k=l  1=1  2h^   ^p 


k 


"5         5(yn   -  X  )5(y     -  X  )  . 
H---is-l=l      Pi  1r 


5(yi    -  X J 


(Montroll's  "interaction  propagator.")  [12] 

^o 
One  of  them,  Fp,  we  already  have  by  summing  (8.14) 

and  (8.16)  : 


oo     ,  ,,    t  t-1   -a(6^-i-a)  + 

(8.20)   F°(k,a)  =^  (t)t+l  ^  A  ^e 


a(p.g4<t) 


t=l 


^=0 


Let  us  now  write  the  partition  function  (8.1)  for  a  system 
in  which  a  t-particle  toron  has  been  fixed: 


(8.21)  ,juj  =:Sf^zz(i)P 


^t-^ 


s-0 


dx,  ...  dx 

1  £ 


dx-,   .  .  .  dx^     exp  -  >   e 

P  ^Pi     ^Pm-1     ^^ 


M     r  t 

^r m   -2 


1 


2^2  ^Ip 


k 


s  s 

,  -v fn   '2   ,1  V —  '  i  /  \ 

27i    Pj^     1        Pj^    Pj^ 


1^' 

2  V 

s        t 

ZIIZ 

My.-  -  y.  ) 
Pk   'Pk 


i=l  j=l 


-K     'K 
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Including  two  constant   terms    (first   and  fourth)    In  y^    ...    y^, 
and  Ui    Is   the   last   term  In  the   bracket.      Then,    according 
to    (8.6)    and    (8.7) : 


00 


t,4.vt+l 


(8.22)      XZ  ^    (-^  /    ^^1    •••    ^^t    ?   ^^1^    '"    ^^t 


t=l 


•Pi 


p. 


M-1 


>[U^]^(yi,,-^'^'^%n-^"^ 


P' 


=  n2(x'P'x"  P"  )p. 


and 


(8.23)     ^  z^    (i)^^-^    /   dy^    •••    dy^j'   cly^      ...    dy^ 


Pn  P 


M-1 


^,[U^]6(y,^^-x')    =n^(.'P')>. 


Similarly  writing  the   density: 


(8.24)      n^    (x"  P"  |U^)^^[U^] 


00     „s 


s=l  p  '-' 


dx.    .  . .    dx^    /      dx-,       .  .  .    dx 

1  =Jp        \  ^Pm-1 


exp   -e 


"  ^     m      -2         ,    T-^     m      -2 


|Zl'  H-.    --i.  )  ^1 


t    , 


1 

s        t 


Pk    -^^k 


<t>(yi    -  y.-    ) 

1q  Jo 

k         ^k 


1  ^P,,  -^P 


ZIXiMx.      -y        ,P    )1         •    6(x  -  X" 


)    , 


k  ^k 
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equation  (8.8)  implies: 

00 

(8.25)  x:-*(t)*"V*^i  •••  ■^^t/^^i,  •••  ^yt« 

\(^"P"|u^)^tfu^]s(yip.-"'^ 

+  ^  z*  (+)*+!  y^dy^  ...  dy^jfdy^   .  .  .  dy 


t=2  J    ^        ^  J    ^6,       ^6 

^1      ^n-1 


t, 

Pi 

t 
^    ^   1=2      6'         ^6"  '^ 


-pi         ipii  2  >-• 

Equations  (8.22)  and  (8.2^)  convert  the  partition  fimction 
to  the  self  pair  distribution  and  density  respectively, 
while  (8.25)  requires  both  partition  function  and  density 
to  convert  to  the  mutual  pair  distribution.   In  any  case 
these  equations  are  the  "bootstrap  operations"  we  set  out 
to  find. 

In  the  next  section,  the  last  devoted  to  path  Integral 
methods,  we  apply  these  results  to  obtain  a  Debye-Huckel 
equation  in  quantum  statistics. 

We  shall  for  convenience  call  the  second  term  in  (8.25) 
n(x'P'x"  P")j,,  that  Is, 

(8.26)  gz*  (±)*^^7>y,  ...    iy^j   dy    ...  dy 

t 

AfU(l)^  ^  ^(^1   -  x')6(y.    -  X"  )  =  n(x'P'x"p"  )i  . 

For  free  particles  (compare  (8.8),  second  term), 

(8.27)  n°  =  n°=  n°(;)  -  n^  . 
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9.   Quantum  Statistics:  Debye-Huckel  Sequence 
Expand  log  j^[Ui]  v.  Ui,  where j^[Ui]  Is  given  by  (8. 19), 


and 


(9-1'   "<>(^i„  'Pk'  =?^*'"lp  -^Jp  '^K^ 


k 


(9.2)   log  ^^[U^]  =  log  ^  JO]  -  J  dP  dx  n^  ^  Mx  -  yj.. 


P) 


0 


4-  i  /   dpdxdpdx 


0 


n  (xP;xP)  +  n  (xPx^)  -  n^ 


— .  I  ,  *• 


J~  H^  -  y.  ,   P)Mx  -  y.^,  P)  + 


Differentiating  (9.2), 


(9.3)   n^  (x"P"  |U^)  =  n-^  -J 


P  __ 
dxdp 


n2(x"P"  ;xp) 


rr   t 


+  np(x"P"xP)  -  n^  y~  M 


,x  - 


P 


Linearizing  (9.2)  and  dropping  the  linear  term  in  (j): 

,P  __ 

^0 


(9.^)  Jt^^^^    =^tf°^  "^  ly   dpdxdpd5^^[0]]^n2(xP;x"p) 


+  n2(xpxp)  -  n^ 


FI  Mx  -  y-i  ,  P)Hx  -  y.^,  P) 
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Applying  (8.22)  and  (8.26),  neglecting  the  t-toron  self 
interaction  on  the  right  hand  side  of  (9.4); 


(9.5)   n2(x'P'x"  P"  )  +  n(x'P'x"  P"  )  =  (n^  +  n)° 

+  |-   /dFd3rdpdxdxdx[n2(xPxP)    +  n  +  n  -  n]^{x   -  x)^{x   -  x) 
0 


•   ^°(x'p'x"  p"xpx^p)    . 


^o 


(See    (8.19)    for  F^. ) 

Multiplying  (9. 3)  ^^y^t'^^.j)^  ^^^  applying  (8.22),  (8.23), 
(8.25)  and  (8.26): 

(9.6)   n(x'P'x"  p"  )  -  n 

p ^  _  __  ^  ^_^ 

dxdpdx[n(x"  P"  xp)  -  n  +  np  +  n](t)(x  -  x) 
0 

' (n2(x'p'xp)+  n)  . 

Equations  (9.5)  and  (9.6)  are  simultaneous,  consistent 
integral  equations  for  n  -  n,  n^  +  n.  The  n  itself  is  to 
a  first  approximation  n 

Just  as  in  Section  7,  we  examine  a  simpler  approximation, 
which  is  (9.6)  with  n  +  n  replaced  by  its  free  particle 
value.   But  then  (9.6)  has  the  same  form  as  (7.8)  with 
kernel  for  the  homogeneous  equation 

-^v  o   t,o 


(9-7)  (Hg  +  n)"  =  Fg 
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This  kernel,  given  by  (8.20),  has  the  properties 

(9.8a)        F°(k,P-a)  =  F°(k,a) 

and 

(9.8b)        F°(k,p-Kx)  =  F^(k,a)  . 

Then  Immediately 

^n   * 
(9.9)      n(k,P'P")  -  n  =  ;^^-2-  ^^O"  )(j,^(p.)  ^ 


n  ^-\ 


with 

(9.10)      ^n  "  /   ^^^  F2^k,a)  exp  2Trln  g  ; 

h     and  A  are  the  same  as  In  Section  7. 

This  Is  again  the  Montroll  and  Ward  result  [10]. 
Equation  (7.I9)  holds  with  the  new  A^.   The  fermlon  solution 
at  low  temperatures  leads  to  the  Gell-Mann-Bruckner  equation 
for  the  correlation  energy  (Coulomb)  of  the  ground  state. 


10.   The  Potential  Ensemble 

The  final  three  sections  of  this  work  are  devoted  to  an 
alternate  approach  to  finding  Integral  equations  for  the 
pair  distribution.   As  we  did  earlier,  we  will  Introduce 
the  main  ideas  by  working  with  the  simpler  formalism  of 
classical  statistical  mechanics,  then  generalize  . 
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Consider  the  classical  grand  partition  function  for 
a  system  with  Interpartlcle  potential  ^   and  external 
potential  U 

00     N     r 

(10.1)  ^[<i)|U]  =  yZ  |r  /  dx-L  ...  dXj. 
^  N=0   *  >-' 

-  ,  N       - 
exp  -  p  ^  J^   (l)(Xj_  -  X  .)  +XZ  U(x.  )   . 

-  i,j  "^    ^~      ^ 

We  will  allocate  the  potential  (j)  to  an  external  potential 
f   which  uses  the  particles  as  sources 

N 

(10.2)  ^(x)  =  YZ   Mx  -  X.) 

j=l         ^ 

(10.3)  =J   dy  Mxy)p(y)  , 

where 

N 

(10.4)  p(x)  =  y~   6(x  -  X  .)  . 

j=l        J 

If  K(xy)  Is  the  matrix  inverse  of  (j), 

(10.5)  /  K(xy)  <|)(yz)dy  =  -  5(x  -  z)  , 

where  the  negative  sign  makes  K  positive  definite  for  attrac- 
tive potentials,  (10. 3)  inverts  to 

(10.6)  /^K(xy)  ^(y)dy  +  p(x)  =  0  . 
The  potential  energy  now  assumes  the  form 
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(10.7)   |ZI  '  ^(^±  -   ^j)  =  /  ^  P(^)  V'(x) 


+  I  /  dxdy  K(xy)^(x)^(y)  , 


apart  from  a  self  term. 

We  anticipate  that  the  partition  function  can  be 
rewritten  as  a  functional  Integral  over  field  amplitudes 
[^(x)j.  I.e. 

(10.8)   Zj^[(i)|U]  =^J   <ix^    ...  dx^JjDfU)    exp(-  p  ^Z  ^(^i) 

N  AT 

-  3  ^~  U(x^))  •  exp(-  I  /  K(xy)^(x)^(y))dxdy 


//  DV'(x)  exp(-  I  /  K(xy)^(x)^(y))dxdy 


when  //  DV'(x)  is  suitably  Interpreted.   The  pure  field 
contribution  existing  In  the  absence  of  particles  has  been 
divided  out. 

That  (10.8)  does  indeed  represent  Z„  can  be  shown  using 
only  properties  of  translatlonal  Invarlance  and  positivlty 
of  the  functional  integrals  [9].   The  meaning  of  these 
integrals  depends  upon  the  choice  of  interparticle  potential 
(and  boundary  conditions)  in  a  particular  problem.   Thus, 
for  example,  a  one  dimensional  Coulomb  gas,  neutralization 
at  origin,  leads  to  a  Feynman  path  integral  [9]. 

Equation  (10.8)  says  that  the  additional  potential  ^ 
is  equivalent  to  an  ensemble  average  over  external  potentials 
^.   In  brief  notation. 
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(10.9)  Zj^[<t>|U]    =   <   Zj^[0|U   +  ^]    >^ 
where 

(10.10)  <  F(^)    >^=JjFiTp)    exp(-  I  <^|K|^>)r^ 


jjexp(-  I  <il^\K\ip>)^^ 

and 

(10.11)  <  ^IkIV'   >  =    /    dxdy  K(xy)    ^(x)^(y)    . 

Equation  (10.8)  Is  for  purely  attractive  ^.      Purely 
repulsive  ^    (K  negative  definite)  is  treated  by  replacing 

f   by  i^. 

The  corresponding  grand  canonical  form  of  (10. 9)  is, 

easily, 

(10.12)  Q  [•t'lU]    =   <^[0|U  +  f]    >^    . 

Differentiating  (10.12)  once  with  respect  to  U: 

(10.15)   n^[x,(|)|U]^[(t)|U]  =  <^[0|U  +  ^]n^[x,0|U  +  ^]  >^ 

(10.1^)   n^[x,<|)|U]  =  <  n-L[x,0|U  +  ^]^[0|U  +  f]    >^ 


<  a[0|U  +  ^]  >^ 

Expanding  the  notation  in  (10.12), 

(10.15)  oA^W    =  JJ   ^0|U  +  f]    exp(-  I  <^|K|^>)D^ 


iT ex-pi-   I  <^|K|V/>)D^ 


JJ        2 
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let  us  find  the  f   which  maximizes  the  Integrand  of  the 
numerator. 

The  "tJJ  must  satisfy 

(10.16)  4=  [2i[0|U  +  f]    exp   -  I  <^|k|^>]    =  0 
or,    using    (10.11) , 

(10.17)  n^[x,0|U  +  iJ]    +  J  dy  K(xy)^(y)    =  0    , 

(10.18)  ^(x)    =    f  dy  n^[y,0|U  +  ^iMxy)    . 

If  we  approximate  n-|^[x,0|u  +  ip]    by  n-|^[x,0lu  +  ^]  In  (10.1^), 
making  it  read 

(10.19)  n^[x,<t)|U]  =  n-L[x,0|U  +  ^], 

the  one  particle  distribution  uncouples  and 

(10.20)  n^[x,0|U  +  f]    =  r\   exp(-  pu(x)  -  PV^(x))  . 

But  Tp     has  been  determined  self-conslstently,  so  that  by 
(10.18), 

(10.21)  log  n^[x,0|U  +  iJ]    =  log  n  -  pu(x) 

-  P  Jdy  n3^[y,o|u  +  ?]<t)(xy)  . 

If  the  external  potential  arises  from  a  particle  fixed  at  z, 

(10.22)  U,(x)  =  <t)(x,z) 
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_   ^  HgCxz) 


(10.23)       n-L[x,0|U^+  ^]  =  -^ 


equation  (10.21)  becomes 


(10.24)   log 


n2(xz)-i 


.  n2  -. 


(xz)  -  I  f  dy  n2(yz)Mxy) 


This  is  the  non-linear  Debye-Huckel  equation  (see  ref.  (1), 
p.  74). 


11.   Lattice  Representation 

Before  the  potential  ensemble  method  can  be  successfully 
applied  to  quantum  systems  It  Is  necessary  to  cast  the 
partition  function  Into  an  appropriate  lattice  form.   We 
begin  with  the  Boltzmann  canonical  partition  fiinctlon 
from  (4.13) : 

(11.1)  Z^  =   A^  /  dx^  ...  dx^  r  dx^   •••  ^N 

^1       ^M-1 

^M    i-N   _   .^     T   N 
exp  - 


M    pN        o     1  _N   I  "I 


The  first  transformation  of  (11.1)  Is  to  make  the  expression 
discrete  In  x  space  as  well  as  In  temperature. 
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Suppose  that  In  each  spatial  dimension  particles  are 
restricted  to  a  grid  of  discrete  locations,  spaced  by  A, 

(11.2)  a^   "  ~  2'  ""  2  "*"  ^''•'>    -   A,0,A,...,  g  '     (each  ^^) 

there  being  Q  lattice  sites  in  all  along  any  line  parallel 
to  a  coordinate  axis,  and 

(11.3)  O  A  =  L  . 
Then 


(11.4)  z,  =  A^NyZ  •••  E 


N 


a,      a,.. 


M 


'M 


m  a2  ^  o 


ij-' 


+  ^JIMa,   -  a   ) 

T7T    Pk   Pk  ' 


(see  4.11).   When  A  — >  0  we  recover  (11.1).   (Normalization 
factors  will  be  omitted  for  simplicity  in  writing). 

One  other  transformation  will  be  made.   Instead  of  the 
position  representation  (11.4),  we  will  specify  the  number 
of  particles  at  a  particular  lattice  site.   Also,  the  lattice 
sites  will  no  longer  exist  in  x-space  alone,  (11.2),  (with  the 
notation  change  a.   — >  x),  but  in  temperature  dimension: 

(11.5)     a  =  p',  p'  -  e,  3'  -  2e,...,  2e,  e,  0, 
M  sites  spaced  by  e. 
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(11.6)  Me  =  P'  ,      (here  P'  =  (kT)"^)  , 


and  particle  number 

(11.7)  1  =  1,  2,.. .,  N  , 

spaced  by  unity.   The  quantity  v(xia)  will  mean  the  number 
of  particles  at  the  site  (x,i,a). 

In  this  "occupation  representation" 

(11.8)  TZ^ia,      -a.  )   =  ^  T~  v(xla)v(yjP)(:x  -  y) 

since  this  interaction  is  iso-particle  number  and  nearest 
neighbor  in  temperature.*  Also 

(11.9)  :y~iTl^Z  '^^A    -^1    )  =  i^ZIZ  w^i°^)v(yj^) 

k^  ^   i   j    ^^^        ^\  xia  yjp 

an  isothermal  interaction. 

Refer  to  the  mixture  analogy  and  diagram  at  the  end  of 
Section  4.   We  do  not  have  to  worry  about  inclusion  of  the 
self  term  on  the  right  of  (II.9).  (1  =  J) >    ^°^   reasons  which 
appear  shortly. 

The  partition  function  transcribes  to 


*  Note:   X  and  a  made  dimensionless,  by  division  by  A,  e 
respectively. 
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(11.10)   z 


N 


J—'  A^  exp  -  e  ^  ^  v(xia)v(yjp) 
■jv(xia)^  xia  yjP 


■"  ^v^^ 


i^7 


(^  -y''«ij5|a-3i,i  +1*'^  -y'^p. 


All  configurations  of  v(xia)  are  summed  subject  to  the 
following  restriction:   for  every  (l,a)  pair  we  must  have 


(11.11) 


v(xia)  =  1  ,  (stronger  than 
^  v(xia)  =  0,  1)  , 


for  particle  i  at  temperature  a  is  to  be  found  at  one  and 
only  one  point  in  space.   The  five  dimensional  grid  we  are 
working  with  plus  the  exclusion  (11.11)  is  shown  in  the 
diagram. 


X 

>  4 

>  ( 

/ 

;i- 

>    ( 

P 

f 

a 


-^  € 


Every  row  of  sites  parallel  to  x  axis  has  one  and  only 
one  occupied  site 


The  total  number  of  occupied  sites  is  MN. 
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This  restriction  to  the  conf Igurational  sum  may  be  dropped 
if  we  append  an  "exclusion  potential"  ^  with  properties: 

(11.12)  f(xyapij)  =  oo    when  i  =  j,  a  =  p  simultaneously, 

=  0     elsewhere. 

We  shall  understand  such  a  potential  to  be  present  In  (11.10). 

The  appropriate  boundary  conditions  on  (11.10)  for  classical 
statistics  are,  for  each  i,  to  fold  the  (x,a)  plane  of  sites 
into  a  cylinder  joined  at  temperatures  (0,P')  so  that  the 
vertical  lines  a  =  0,  a  =  P'  share  the  same  sites.   One  of 
the  virtues  of  the  method  to  be  developed  is  that  boundary 
conditions  need  never  be  taken  into  account. 

We  next  generate  distributions  in  the  usual  way: 
append  an  external  potential  and  differentiate. 

Taking 

(11.13)  U  =  _2^  v(xia)U(xia)  , 

xia 

5  log  Z  [U] 

(11.14)         tl =  n,(x'i'a' |U)  , 

-6U(x'i'a')      ^ 

bn, (x'i'a" |U) 

(11.15)  — =  no(x'i'a',xla|U) 

-6U(x  i  a)       "^ 


+  n-,  (x'i'a' |U)6  ...  ,  -7-  a""*" 
1^       '  '  x'i'a', xia 

-  n-^(x'i'a' |U)n^(xia|U)  , 


where 


(11.16)  n  (x'l'aM  =  Zn\?II  .,  ^^"^  v(x'i'a')  exp  -  [  ]  , 

1  \^v{xla  )3 

and 

(11.17)  n,(x'i'a',Ha)  =  ^~^  TZ       A^'^  v(x' I'a' )  v(xla)exp-  [  ]  . 

^  iv(xla)j 

These  functions  are  related  to  the  space-temperature  distribu- 
tions of  the  path  integral  methods:  taking  lim  , 

N 

(11.18)  5   n3^(xla)  =  n^(xa)  , 

(11.19)  ^  [n2(xia,yip)+  n^i^^a)b^^^^^'   A"^]  =  n^i^ay^) ,    "  self", 

(11.20)  ^^n2(xia,yJP)  =  n2(xa;yP),   "mutual"   . 
i=l  j=l 

The  second  term  in  (11.19)  Picks  up  the  singularity  in  the 
self  distribution  when  a  =  P,  (see  discussion  after  (3.I3)). 

Bootstrap  Operation 

The  density  in  an  external  potential  due  to  a  single 
occupied  lattice  site  is  proportional  to  the  pair  distribu- 
tion for  a  uniform  system.   To  see  this  Isolate  an  occupied 
site  at  xia,  v(xI5)  =  1.   This  site  supplies  an  external 
potential  to  the  others: 
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(11.21)    Uj  j^  =  y^   v(xla)v(xla) 


xia 


r   2 


m^  U    -D'   5i^i  5|a-5|,l  +  ^(^  -  ^)5^,a 


The  density  in  the  presence  of  Ui_^j^  is: 

I        'NIIVI  — 1  ' 

(11.22)      n    (x'l'a'IU,    j^)    ==  IZ  ,    A      "      v(x'l'a')exp   -    [    ]    exp-eUi^^^ 
-^  ^  ^v(xia)\ 


'       NM  ' 

TI       A        exp   -    [    ]      exp-eU,    ^ 

^v(xia )] 

the  prime  on  the  exponential  meaning  that  the  site  v(xla)  Is 

excluded  from  the  sums  appearing  there.   This  site  enters  Into 

Ui  ^  only.   Similarly,  the  prime  on  the  configuration  sum 

means  the  sum  is  taken  over  a  lattice  in  which  the  one  site 

at  (xia)  has  been  excluded.   Dividing  top  and  bottom  by  A 

in  (11.22),  we  see  that  this  is  precisely  the  meaning  of 

n  fx'i^ct'  xiot) 

-  defined  on  the  full  lattice  -  i.e.  augmented  now 


nn (xia) 

by  (xIa).   So 

np(x'i'a' ,xia) 

(11.25)      n  (x'i'a'  I  U,  ^)  =  -— • 

J-  y-Ti^  n.(xia) 

This  relation  has  the  simplicity  of  the  classical  one,  see  (2.l6) 
It  is  just  at  this  point  that  the  relation  to  the  path  integral 
methods  becomes  apparent:   had  we  fixed  an  array  of  occupied 
sites  one  per  a  value,  a  =  0, . . . ,  M,  (1  =  constant),  the 
resulting  bootstrap  equations  would  be  the  same  as  in  Section  5. 
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The  next  step  has  been  to  expand  directly  n-|^[U]  v.  U  in  a 

Taylor  series.   Here  we  are  in  a' dilemma.   If  we  expand  in 

-1       2 
the  full  Ui^j^g,  the  KE  is  so  strong  (-'e   (x  -  y)  )  it  over- 
powers the  pair  interaction  resulting  in  nonsense,  while  if 
we  expand  in  U,,  it  is  too  weak,  (--'e),  to  give  any  correction 
to  free  particle  solutions.   The  potential  ensemble  removes 
these  difficulties  by  replacing  the  pair  potential  by  an 
external  potential  ^  determined  self- consistently. 


12.   Application  of  Potential  Ensemble  to  the 
Lattice  Representation 

Working  with  a  more  general  pair  potential  <j)(xla,yjP) , 
and  with  kinetic  energy  and  exclusion  potential, 

2 

(12.1)  0(xia,yjP)  =  ri^  (x  -  y)^  5ij  ^ja-pUl  +$(^la»yJP)  ■ 

(11.10)  becomes,  in  brief  notation: 

(12.2)  z^[e   +  ^|U] 

=  y-     A^  expf-  e  XZ  v{x)v(y)[i.  0(xy)  +  ^  H^y)]] 


exp  -  e  >   v(x)U(x)  . 

X 
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We  have  added  an  external  potential,  and  used  the 
notation  x   =  (xla),  y  =  (yjP). 

Allocate  the  ^   to  an  external  potential  ^  using 
occupied  lattice  sites  as  sources: 


(12.3) 


^(x)  =  y^  <[>(xy)v(y) 
y 


If  K(xy)  is  the  matrix  inverse  of  H^y) 


(12.4) 


so  that 


">       K(xy)(t)(yz)  =  -  5 


y 


xz  ' 


(12.5) 


5   K(xy)^(y)  +  v(x)  =  0  , 

y 


the  potential  energy  assumes  the  form 


(12.6)  |XI  v(^)^(y)Mxy)    =  y~  v(x)V/(x) 


x,y 


+  |lZ  K(xy)  ^(x)^(y)    . 
x,y 


In  analogy  to   (10.8), 


(12.7)      Z^[e   +  <t)|U]    =;^     A^^   /Texp    -    e  2Z  v(x)v(y) 

(v { X )]  -^  x,y 


e(^y) 

2 


exp    -   e   >       v(x)(V/(x)    +  U(x)) 

X 


exp   -(#  3~  K(xy)^(x)^(y))D^(x) 


x,y 


/7exp(-  I  XI  K(xy)V(x)^(y))D^(x) 

'-iw  X  .V 


x,y 
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briefly. 


(12.8)         z^[e  +   ^|U]  =  <  Zj^[0|U  +  f]    >^ 


A  differentiation  with  respect  to  U  yields: 

(12.9)   n^(x,e  +  (|)|U)  =  <  n^{x,e|U  +  ^)Z^[0|U  +  f]    >^ 


<  z^[e|u  +  ^]  >^ 

The  -ip   which  maximizes  the  numerator  of  (12.8)  is  given  by 

(12.10)        ^(x)  =  2Z  ^<i>(^y)ni(y^^|u  +  ?)  , 

y 

V 

compare  (10.18) . 

Now  let  us  approximate  n-,  (x,elu  +  -ip)    by  n-,(x,0|U  +  ^) 
in  (12.9)  making  it  read 

(12.11)  n^(x,0  +  ({)|U)  =  n^(x,0lU  +  ^)  . 

Suppose  now  that  U  arises  from  an  occupied  lattice  site  at  z, 
so  that  it  is  composed  of  two  parts,  kinetic  and  potential 
energies : 

(12.12)  U  =  Uq  +  U^  . 

Expand  n, (x, 0|Ug  +  Ui  +  ^)  v.  Ui  +  ^,  using  (II.I5), 
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(12.13)  n^(x,e|Ug  +  \J^   +  f)    =   n^(x,e|U0) 

-  2Z  [np(xy,  e|Ug) 

y 

+  n(x,elUg)  ^^ 

-  n(x,0lU0)n(y,0|U0)]eA 

(U^  +  V^){y)  +  ...   . 

But  ^(y)  has  been  determined  self- consistently  by  (12.10), 
so  substituting  this  into  (12.15), 

(12.14)  n^(x,0|U0  +  U^  +  ^) 

=  n  (x,e|U  )  -  XI  [n  (xy,0|Ug)  +  n(x,0|U0)  -^ 

y 

-  n(x,0|Ug)n(y,0|U0)]eA 
•  [Ui  (y)  +  ^  A<l)(ycD)n^(a),0|Ug  +  U i  +  ^)  ] 

Now  Ui  +  ^  is  tied  to  the  fixed  occupied  site  at  z, 

(12.15)  u^(y)  =  Hyz)  , 

and 

(12.16)  n^{^,e\\jQ   +  U^  +  ?)  =  ^^§^j^ 

using  (11.23)  and  (12.11). 

Equation  (11.25)  also  means  that 
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Or  \ 

n  (  x  z  ) 
(12.17)       n^(x,e|Ug)  =  2 


n°(z) 


and  by  a  simple  generalization  of  (11.25), 

n^(xyz) 

(12.18)       n  (xy,  e|Ug)  =  -2-- . 

n-|^(z) 


Substituting  all  of  these  into  (12.14), 

[n^(xyz)  +  n2(xz) 


(12.19)   n^(xz)  =  n°(xz)  -  ^Z   [n°(xyz)  +  n°(xz)  ^^ 


Of         \ 

o        ^o(y^^        1 

•  eA^  (t)(yco)[n2(cDz)  +  n^(z)  -|^  ] 


f 


an  integral  equation  for  np(xz).   One  equation  suffices,  for 

if  n_(xz)  is  found,  the  self  and  mutual  pair  distributions 

with  space  and  temperature  arguments  only  are  gotten  via 

(11.17)  -  (11.19). 

If  we  perform  a  superposition  approximation  in  the  first 

bracket  of  (12.19),  separating  out  the  z  dependence: 


5 
(12.20)   np(xz)  =  n°(xz)  -  ^  [n°(xy)  +  n  (x)  -|^  -  n, (x)n, (y) ] eA' 

•<t)(yo^)[n2(coz)  +  n-^(z)  -^]  . 

The  potential  is  isothermal  and  has  no  dependence  on  particle 
number.   So  that  expanding  the  notation: 


7^^ 


(12.21)  n„(xia,zkv)  =  n°  -  XZ  ^  [n°(xla,yjp) 

^  "^   y,a3  ^~ 

+  n(xla)5^^^^yjp  A"   -  n(xla)n(yjP) ] 
•  EA^  (t)(y  -  cd)  [n2(cDip,zkv)  +  n(zkv)A"-'- 

^ai£p,zkv^  • 

The  easiest  way  to  compare  this  to  earlier  work  is  to  sum 
1,  j,    k  and  &    independently  from  1  to  N.   Then  according  to 
(11.17)  -  (11.19),  e,  A  — >  0, 

(12.22)  n^{:>^a;zv)    +n2(xazv)  =n°{;)    +  n°{    ) 

-    /  dydcDdp[n°(xa;yp)    -  n^   +  n°(    )]\){y   -  co) 

•    [n^iu^^;zv)    +  n^{    )  ]    . 

Replacing  np(  )  by  n°(  ),  leaves  us  with  (7-8)  in  classical 
statistics  or  (9-6)  in  quantum  statistics  (with  n' replaced  by 
n°  in  (9.6)). 

The  method  developed  here  puts  both  statistics  on  an 
equal  footing  because  the  boundary  conditions  on  Z^  and  the 
distributions  never  enter. 

It  is  not  really  necessary  to  make  the  separate 
approximation  n^{    )    =  n^{    ).   If  (12.22)  is  solved  this 
approximation  will  result  when  the  thermodynamic  limit 
(N  — >  oo,  V  — ^>  OD  ,  TT  =  n)  is  taken. 

Equation  (12.19)  is  an  improved  approximation  over 
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(12.20).   Also,  we  could  carry  out  the  Taylor  expansion 
(12.13)  to  higher  orders.   In  addition  one  can  improve  the 
steepest  descent  calculation  represented  by  (12.11).   To  do 
this  go  back  to  (12. 9)  and  expand  n^(x,0|U  +  ^)  about  f: 

(12.25)   n.(x,e|U  +  ^)  =  n  (x,0lU  +  f)    -   ^Z  ^^2^^^'   ^'^  ■*"  "^^ 

y 

+  n^(x,e|U  +  f)    -^ 

-n-^(x,0|U  +  lp)n^{y,e\u  +  ^)] 

'  Ae{ip{y)    -  f{y) )  +  . . . 

The  rir^   in  the  bracket  can  be  expressed  in  terms  of 
n, (x,0lu  +  ^)  by  a  Taylor  expansion  similar  to  (12.13)  with 
the  self-consistency  relation  (12.10)  employed  once  more. 
Expectations  over  powers  of  tp   would  have  to  be  done  (requiring 
a  specific  <j)  and  boundary  conditions)  to  improve  upon  (12.11). 
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Appendix 
Notation   For    Path   Integral   Methods 


n   ,    distributions 
s 

A 

n  ,  modified  distributions 

s ' 

F  ,  Ursell  Functions 
s' 

K 

F  ,  modified  Ursell  Pimctions 


np(  ),  np,  self  pair  distribution 
np(;),     mutual  pair  distribution 


n  =  n^ ( ; )  -  n 


n-,  ,  n,  density 
n,      see  (8.26) 


°  =  "  free  " 

Hplk,  a),  Fourier  transform,  self  pair  distribution 
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